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ORIGINAL ARTICLE 

Epigenetic modulation of glucocorticoid receptors in 
posttraumatic stress disorder 

B Labonte 1,2,5 , N Azoulay 2 ' 5 , V Yerko 1 , G Turecki 1 ' 2 ' 3 ' 4 and A Brunet 2 ' 3 ' 4 

Some individuals suffering from posttraumatic stress disorder (PTSD) exhibit lower basal salivary Cortisol and higher glucocorticoid 
receptor (GR) sensitivity. Recent studies suggest that epigenetic mechanisms regulate the activity of Cortisol and GR. As a means to 
combine and cross-validate those findings, we compared Cortisol, GR expression and promoter methylation levels in peripheral T 
lymphocytes of healthy controls versus individuals endorsing a diagnosis of lifetime PTSD. Thirty subjects with lifetime (current or 
remitted) PTSD and 16 subjects never exposed to trauma were recruited. Salivary Cortisol was collected at six time points over the 
course of a single weekday and analyzed utilizing a time-resolved fluorescence immunoassay. GR expression (GR to tai/ 1b# 1o 1f and 
1 H ) was measured by quantitative RT-PCR. DNA methylation levels in human glucocorticoid receptor (hGR) 1 B and 1 C variant's 
promoter were quantified by epityper in T lymphocytes isolated by magnetic-assisted cell sorting. Individuals with lifetime PTSD 
have lower morning Cortisol release, higher mRNA expression of hGR to tai/ 1b/ and 1 C and lower overall methylation levels in hGR 1 B 
and 1 C promoters. Cortisol levels were inversely correlated with hGR 1 B mRNA expression. Moreover, overall and CpG site-specific 
methylation levels were inversely correlated with hGR tota i and 1 B mRNA expression. There was no difference between current and 
remitted PTSD across Cortisol, GR expression mRNA and DNA methylation data. Traumatic events induce DNA methylation 
alterations in distinct promoters of hGR with transcriptional modifications that associate with hypoactive hypothalamus-pituitary- 
adrenal axis in individuals with PTSD. Our results also point toward an important role of hGR 1 B variant in PTSD. 
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INTRODUCTION 

Posttraumatic stress disorder (PTSD) is characterized by symptoms 
of intrusions, avoidance and hyperarousal following exposure to a 
traumatic event. Lifetime trauma exposure has been estimated 
between 40-90%. 1-6 However, the lifetime prevalence for PTSD is 

7-1 2%y 

If remission is considered the normative response to trauma 
exposure, a substantial proportion of individuals go on developing 
long-lasting intractable PTSD. The clinical manifestations of PTSD 
may reflect enduring changes in the stress-related neurobiological 
structures and systems. 8 Indeed, changes in the activity of the 
hypothalamus-pituitary-adrenal (HPA) axis, the main stress reg- 
ulatory system, are induced by exposure to severe stressors in 
animals and humans 9,10 and dysfunction of the HPA axis activity is 
a well-characterized feature in PTSD. For instance, such patients 
exhibit higher corticotropin releasing factor in the cerebrospinal 
fluid, 11,12 lower plasma and serum morning Cortisol release, 13 
higher glucocorticoid receptor (GR) expression in lymphocytes 14 
and greater GR sensitivity. 15 In addition, recent studies have 
shown that life-events, at least early in life, can modify the 
chemical structure surrounding the genome and influence gene 
expression, 16 including GR expression. 17-19 

The human GR gene contains eight coding exons (exons 2-9) 
and nine 5' non-coding first exons (exons A, B, C, D, E, F, H, I, J) 
spanned through a region of more than 80 kb within chromosome 



5q31-q32. Both rodent and human GR levels are regulated 
almost exclusively at the transcriptional level 24,25 and each non- 
coding exon 1 variant has its own promoter 22,26 " 30 These multiple 
independent promoter regions contribute to the tissue specificity 
of GR expression 23 DNA methylation is one of the epigenetic 
mechanisms controlling the expression of genes and has been 
shown to be modified by environmental stress. 31 For instance, in 
rats, maternal behavior regulates hippocampal GR1 7 expression in 
offspring through changes in GR1 7 promoter methylation. 17 
Similarly, decreased hippocampal total GR, 1 B , 1 C/ 1 F (rat homolog 
for human GR1 F ) and 1 H expression associates with respective GR 
exon 1 variant promoter methylation alterations in suicide 
completers who were abused in childhood. 18,19 These results 
suggest that methylation changes in GR promoters modify GR 
expression in individuals who experienced early-life adversity. It 
remains unknown if epigenetic regulation also explains HPA axis 
activity in other conditions associated with adversity such as PTSD. 
However, by definition, PTSD is etiologically related to environ- 
mental stress and, as such, it is likely to involve an important 
epigenetic component. 

To address this question, we examined HPA axis functioning in 
individuals endorsing a diagnosis of lifetime PTSD by investigating 
salivary Cortisol measures, the expression of four variants of GR (1 B , 
1 C , 1 F and 1 H ), and their association with variations in exon 1 
variant promoter DNA methylation in T lymphocytes. We 
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investigated T lymphocytes because these cells have an important 
role in cell-mediated immunity, and as such, are depositaries of 
cellular memory to environmental factors. In addition, DNA 
methylation patterns vary between different cell fractions, 32 thus 
it is preferable to investigate DNA extracted from individual cell 
populations rather than from homogenates. We hypothesized that 
the HPA axis alterations in PTSD are associated with variations in 
Cortisol levels due to changes in the expression of the GR gene. 
We further hypothesized that modifications in DNA methylation 
levels within the promoter region of the GR gene are associated 
with changes in GR expression in PTSD. 



MATERIALS AND METHODS 

Subjects and procedure 

Twenty seven individuals meeting DSM-IV criteria for lifetime PTSD 
(Current PTSD: N = 15; Past PTSD: N = 12) were recruited. Diagnosis was 
established with a semi-structured interview 33 and further corroborated by 
a self-report symptom measure 34 Sixteen individuals with a negative 
trauma history were also recruited. Current comorbid depressive 
symptoms were assessed with a subscale of the symptom check list-90 
using a cutoff of 70 for clinically relevant symptoms of depression. 35 The 
index traumatic events were as follows: motor vehicle accident (n = 19), 
participation in a peacekeeping mission (n = 1), assault with a weapon 
(n = 2), physical and/or sexual abuse in adulthood (n = 14), physical and/or 
sexual abuse during childhood (n = 13) or other (n = 8). All events met the 
DSM-IV Al (life threat) and A2 (fear, helplessness or horror) criteria for 
trauma exposure. Exclusion criteria for controls were a lifetime trauma 
exposure, a current DSM-IV axis I psychiatric diagnosis 36 including any 
form of substance abuse and living outside of the greater Montreal 
metropolitan area. Table 1 lists sociodemographic and clinical sample's 
information. The Research Ethics Board at the Douglas University Mental 
Health Institute approved the project. Signed informed consent was 
obtained from every participant. 



Salivary Cortisol 

Basal salivary Cortisol (unbound free Cortisol) was sampled following 
previously described methods. 37,38 Participants were given six saliva 
collection tubes (Sarstedt, Montreal, QC, Canada ) to use at home over the 
course of a single weekday at the following times: awakening, awakening 
+30min, awake ning+45 min, awakening+60 min, awakening+8 h and 
bedtime. The saliva samples were stored at - 80 °C. Samples were analyzed 
in duplicate all at once using a time-resolved fluorescence immunoassay 38 
The average coefficient of variation between duplicates was 5.2%. Cortisol 
data were plotted on a graph and the area under the curve was measured 
for each time interval according to previously described methods 39 



Peripheral blood RNA extraction 

Blood samples drawn for the RNA extraction (12 ml) by a registered nurse 
were immediately put on Leucolock (total RNA isolation system, Ambion, 
CA, USA) and frozen at -80°C until RNA extraction the next day. As 
opposed to DNA (see next section), RNA was extracted from whole blood. 
Once extracted, concentrations and quality were measured using a 2100 
BioAnalyzer (Agilent Technologies, Palo Alto, CA, USA). Samples with a RNA 
Integrity Number greater than 5.5 on a scale of 1 (lowest) to 10 (highest) 
were used for further analysis. 

Human glucocorticoid receptor expression analysis 
RNA was converted to cDNA using oligo dT primers according to 
manufacturer's recommendations (Invitrogen, Burlington, ON, Canada) 
and expression was analyzed using a real time polymerase chain reaction 
(RT-PCR) in an ABI 7900HT (Applied Biosystems, Foster City, CA, USA). 
Amplification of cDNA was performed using custom-designed taqman 
probes used in a previous study (Supplementary Table SI). Primers and 
taqman probes were obtained from Integrated DNA Technology (San 
Diego, CA, USA; Supplementary Table S1). Quantification was performed 
with five replicates each in a 12ul reaction as followed: 0.6 ul of GR 
primers/probes, 3.1 ul of water, 6.0 ul of Taqman Master Mix and 2 ul of 
cDNA (1 ng ul -1 ). Thermal cycling conditions were: 50 °C for 2 min, 95 °C 
for 10 min and 40 cycles, each of 95 °C for 15 s and specific melting 
temperature for 60s (Supplementary Table S1). The cycle threshold values 
for all replicates were pooled to obtain a mean value. Replicates with a 
standard deviation greater than 0.3 were excluded from the analysis. For 
each replicate, a quantity was extrapolated from a standard curve 
including six dots, each of which was five times diluted (5:1). A standard 
curve was composed of a mix cDNA from all subjects. Mean quantities 
from all sample replicates were normalized to each of the reference genes 
((3-Actin, GAPDH and 18S) and averaged together. Expression data are 
expressed as ratios of average mean quantities. 

T-lymphocyte isolation, DNA extraction and methylation 
quantification 

Blood samples drawn for the T-lymphocyte isolation and DNA extraction 
(12 ml) were immediately put on ice before being processed on the same 
day. T lymphocytes were isolated following the Automacs Pro protocol 
using CD3 whole blood microbeads according to manufacturer's 
recommendations (Miltenyi Biotech, Boston, MA, USA). T-lymphocyte 
DNA was extracted using the DNA micro kit (Qiagen, Toronto, ON, 
Canada). DNA was diluted to 500 ng and sodium bisulfite conversion was 
performed using the EpiTect Bisulfite kit (Qiagen) according to the 
manufacturer's instructions. Bisulfite-treated DNA was then sent to the 
Innovation Center of Genome Quebec (Montreal, Quebec, Canada) where 
epityper (Sequenom, San Diego, CA, USA) was performed on the regions of 
human glucocorticoid receptor (hGR) promoters defined in Supplementary 
Figure S3. Methylation data are expressed as mean percentage of overall 
methylation (global methylation) and methylation at each CpG site. 

Statistical analyses 

Analyses were conducted using SPSS 19.0 (IBM Corp., Armonk, NY, USA). 
Analyses were adjusted for age and sex. Gene expression values were 
normalized and compared using univariate analysis of variance (ANOVA). 
Cortisol and methylation levels were analyzed using a mixed model analysis 
with group (PTSD and controls) as a fixed factor and Cortisol sampling time 
(six samples) or CpG sites as a repeated measure, followed by LSD post hoc 
tests. For each analysis, outliers within groups were excluded when their 
value differed by more than two SDs from the group mean (Z-test). 
Correlations between Cortisol, expression and methylation levels were 
analyzed with Pearson's correlation. Statistical significance was fixed at 
P < 0.05. Data are presented as mean ± standard error of the mean (s.e.m.). 

RESULTS 

Salivary basal Cortisol levels 

A mixed two-way ANOVA performed on overall Cortisol levels 
(area under the curve) revealed a significant main effect of group 
(F(i,87) = 15.2, P< 0.0005), with Cortisol levels being lower in the 
PTSD group compared with controls (Figure 1). The analysis also 
revealed a main effect of Cortisol sampling time (F (4/4 o) = 151.1, 



Table 1. Demographic and samples' clinical information 




Lifetime PTSD 


Controls 


Male:Female 


15:15 


8:8 


Age (s.d.) 


43.4 (12.1) 


36.9 (11.2) 


Type of trauma 






Motor vehicle accident 


19 




War 


1 




Armed assault 


2 




Adulthood physical/sexual abuse 


14 




Childhood physical/sexual abuse 


13 




Other 


8 




Probable depression 


36% 


0% 


Current psychotropic medication use 


40% 


0% 


RNA integrity number (s.d.) 


7.1 (1.3) 


7.2 (0.6) 


Abbreviation: PTSD, posttraumatic stress disorder. PTSD symptom score is 
based on the self-report IES-R that ranges from 0 (no symptom) to 88. 
Probable depression was determined by a score >70 on a subscale of the 
SCL-90. 
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P< 0.0001) and a group by time interaction (F (4j4 o) = 4.6, 
P< 0.005). Post hoc tests indicated that Cortisol levels were signi- 
ficantly lower in the PTSD group compared with the controls for 
the time interval between 30 and 45 min post awakening (P< 
0.01) and between 45-60 min post awakening (P< 0.05; Figure 1). 
In addition, a between-groups trend was found for the time 
interval between awakening and 30 min post awakening (P = 
0.058). An analysis of covariance (ANCOVA) revealed no between- 
group effect of age, sex and depression on Cortisol levels. 

We also compared Cortisol levels in current and past PTSD 
patients. Interestingly, we found no significant difference in the 
Cortisol levels between current and past PTSD (F (4;2 3) = 0.1 06, 
P = 0.98) suggesting that changes in the Cortisol levels are main- 
tained after remission from PTSD symptoms (Supplementary 
Figure S1 A). These effects were also independent of age and sex. 



Whole blood human glucocorticoid receptor expression levels 
The relative expression of hGR to tai as well as the first non-coding 
exons 1 b# 1 o 1 f and 1 h was compared across the PTSD and control 
groups in peripheral blood samples. Results showed a significantly 
higher expression of hGR tot ai (F (1/43) = 5.3; P<0.05; Figure 2a), 1 B 
(F (l 43) = 12.1; P< 0.005; Figure 2b) and 1 c (F (1 , 42) = 4.2; P<0.05; 
Figure 2c) in the PTSD group compared with controls. In contrast, 



Salivary Cortisol Release 



16 




PTSD 
CTRL 



Wake up +30mins +45mins +60mins 4:00 PM Bedtime 

Figure 1. Diurnal variations in basal salivary Cortisol levels (mean ± s. 
e.m.) at awakening, +30 min, +45 min, +60 min, 4pm and bedtime as 
a function of group status. Note that time intervals are plotted 
evenly. Lifetime PTSD: n = 22, Controls: n = 15. *P<0.05; + P = 0.058. 
CTRL, control; PTSD, posttraumatic stress disorder. 



hGR 1 H relative expression was not significantly different between 
the groups (Figure 2d). In addition, consistent with previous 
findings showing that hGR 1 F expression is mainly restricted to the 
brain, 23 we were unable to quantify its expression in peripheral 
blood. ANCOVA revealed no effect of age or sex on hGR tota |, 1 B 
and 1 c expression levels between groups although we found a 
significant effect of age on hGR 1 H expression. 

Interestingly, we found a significant negative correlation 
between hGR 1 B relative expression and Cortisol levels for the 
time intervals between awakening and 30 min post awakening 
(r = -0.34, P<0.05; Supplementary Figure S2A), 30-45 min post 
awakening (r = -0.39, P<0.05; Figure 3a) and 45-60 min post 
awakening (r = -0.37, P<0.05 Supplementary Figure S2B). There 
was no correlation between hGR tot ai or hGR 1 c relative expression 
and Cortisol levels. 

We next assessed whether hGR non-coding exons' expression 
differs in current and past PTSD individuals. Our analyses revealed 
no significant difference in the expression of hGR tota i (F 
( 1 , 26) = 2.277, P = 0.145), 1 B (F ( i,26)=1.611, P = 0.217), 1 C (F (lf25) = 
0.102, P=0.752), and 1 H (F (1#2e) = 1 .1 71, P = 0.290) between current 
and past PTSD (Supplementary Figure S1B). These effects were 
also independent of age and sex. 

T-lymphocyte methylation levels 

We sought to determine whether alterations in DNA methylation 
patterns within hGR variants' promoter could explain the 
expression changes observed between the PTSD group and 
controls. We focused on the two variants (hGR 1 B and 1 C ) that 
showed higher expression in the PTSD group compared with 
controls (Supplementary Figure S3). 

Human glucocorticoid receptor 1 B 

Methylation levels were analyzed across a region of 479 bp 
composed of 29 CpGs located upstream hGR 1 B (Supplementary 
Figure S3). Mixed two-way ANOVA revealed a significant main 
effect of group on percentage of methylation (F (lfl82 ) = 28.7, 
P< 0.0001), CpG site (F (15;96) = 307.8, P< 0.0001), and a group by 
CpG site interaction (F (15 96) = 2.8, P < 0.005). Our analysis indicated 
a significant effect of age (F (1 , 52 ) = 9.5, P< 0.005) on methylation 
levels although the main effects of group, CpG site, and group by 
CpG site interaction remained significant after correction for age. 
In accordance with the expression data, overall methylation level 
in hGR 1 B promoter was significantly lower in PTSD compared with 
control (Figure 4a). Posf hoc tests showed site-specific 
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Figure 2. Glucocorticoid receptor mRNA expression levels (average relative mean quantity ± s.e.m.) of (a) hGR tota |, (b) hGR 1 B , (c) hGR 1 c , and 
(d) hGR 1 H as a function of group status. PTSD: n = 27, Controls: n = 16. *P< 0.05. CTRL, control; hGR, human glucocorticoid receptor; PTSD, 
posttraumatic stress disorder. 
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Figure 3. (a) Correlation between hGR 1 B mRNA relative expression and salivary Cortisol levels for the time interval between 30-45 min after 
awakening, (b) Correlations between total DNA methylation level in GR1 B promoter and hGR 1 B (gray dots) and hGR tota | mRNA relative 
expression (black dots), (c) Correlations between hGR 1 B mRNA relative expression and DNA methylation level in hGR1B promoter at CpG 
sitesl 1 (gray dots) and 13 (black dots), (d) Correlation between Cortisol levels for the time interval between 30-45 min after awakening and 
hGR1B promoter DNA methylation levels at CpG sites 11 (gray dots) and 13 (black dots). hGR, human glucocorticoid receptor. 
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Figure 4. Overall percentage of methylation (mean ± s.e.m.) in (a) 
hGR 1 B/ and (b) hGR 1 c promoters as a function of group status. 
PTSD: n = 30, Controls: n= 16. *P< 0.05. CTRL, control; hGR, human 
glucocorticoid receptor; PTSD, posttraumatic stress disorder. 

hypomethylation at CpG sites 2-4 (P < 0.005), 1 1 (P < 0.05), 1 3 
(P<0.01) and 14 (P = 0.05) in individuals suffering from lifetime 
PTSD compared with controls (Supplementary Figure S4A). 

Our results show a negative correlation between overall 
methylation levels in the hGR 1 B promoter and hGR tota | (r = -0.35, 
P<0.05) and 1 B (r=-0.51, P< 0.001) relative expression 
(Figure 3b). GR to tai and 1 B relative expression was also negatively 
correlated with methylation levels at CpG sites 2-4 (hGR to tai- 
r = -0.37, P<0.05; hGR 1 B : r = -0.55, P< 0.005; Supplementary 
Figure S2C), 11 (hGR 1 B : r = -0.34, P<0.05), and 13 (hGR 1 B : 
r = -0.32, P<0.05; Figure 3c). In addition, methylation levels at 
CpG sites 11 and 13 were positively correlated with Cortisol levels 
obtained between awakening and 30 min post awakening (CpG 
site 11: r=0.4, P = 0.05; CpG site 13: r = 0.44, P<0.01; 
Supplementary Figure S2D), and +30 min and +45 min (CpG site 
11: r = 0.4, P<0.05; CpG site 13: r = 0.43, P<0.01; Figure 3d). 

We assessed DNA methylation levels in current and past PTSD 
samples and found no significant group difference (Supplementary 
Figure S1C). These effects were also independent of age and sex. 



Human glucocorticoid receptor 1 c 

Methylation levels were analyzed across a 364 bp region 
composed of 54 CpG upstream hGR 1 c (Supplementary Figure 
S3). Mixed two-way ANOVA revealed a significant main effect of 
group on percentage of methylation (F (1 ,3 54 ) = 6.5, P<0.05), CpG 
site (F (21 , 69 8) = 245.6, P< 0.0001), but no group by CpG site 
interaction. In accordance with the expression data, total 
methylation level in hGR 1 c promoter was significantly lower in 
PTSD compared with the controls (Figure 4b). Posf hoc analysis 
showed significant hypomethylation in PTSD patients at CpG sites 
40-41 (P<0.05). We also found one CpG site (CpG 51) to be 
significantly hypermethylated (P<0.05; Supplementary Figure 
S4B). An ANCOVA revealed no effect of age and sex on hGR 1 c 
methylation levels between groups. 

Our analyses showed no correlation between total methylation 
in 1 c promoter and relative hGR 1 c expression or Cortisol levels. 

We also compared DNA methylation levels in current and past 
PTSD samples. Our analysis revealed no difference in DNA 
methylation levels between current and past PTSD (Supplemen- 
tary Figure S1D). These effects were also independent of age 
and sex. 



DISCUSSION 

In this study, we investigated the HPA axis profile of individuals 
suffering from lifetime PTSD and healthy controls by measuring 
basal salivary Cortisol levels, hGR peripheral expression and CpG 
methylation states in T lymphocytes. We focused our investigation 
on the non-coding exons 1 B , 1 c , 1f and 1 H , which have been 
previously studied in the brain of depressed subjects 40 and 
abused suicide completers. 18,19 By combining these three 
molecular measures, we sought to obtain a comprehensive view 
of the HPA axis profile and to investigate whether epigenetic 
mechanisms, in response to traumatic events, modify GR 
expression and Cortisol release in individuals with a diagnosis of 
lifetime PTSD. 
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In agreement with the previous findings, 1 3,1 5,41,42 our results 
suggest a dysregulation of the HPA axis activity in individuals with 
PTSD that was independent of current versus past PTSD status. 
Accordingly, the rise in morning Cortisol release, up to 60 min after 
awakening observed in the control group was lower in the PTSD 
group. Importantly, the measure of Cortisol in this study was 
obtained from saliva, which has been highly correlated with serum 
unbound Cortisol levels. 43 Under normal conditions, unbound 
Cortisol increases by about 50% within the first 30 min after 
awakening. 44 Once released in the blood, 90-95% of Cortisol binds 
corticosteroid binding globulin proteins whereas the remaining 
fraction enters the cells and activate both the glucocorticoid and 
mineralocorticoid receptors. 45,46 The lower levels of morning 
salivary Cortisol found in this study may indicate that the HPA axis 
in individuals with lifetime PTSD may be hypoactive. One may 
assume that this hypoactivity is the consequence of GR 
hypersensitivity. 47 

Accordingly, our data showed that hGR expression levels are 
increased in lifetime PTSD independently of current PTSD status. 
Although similar results have previously been reported, 48 it is the 
first time that individual hGR variants are investigated in PTSD, 
and our results suggest that the higher hGR expression may be 
attributed to the cumulative elevations in hGR 1 B and 1 c 
expression. On the other hand, hGR 1 H expression was stable 
between groups and hGR 1 F could not be detected. This is in 
agreement with previous findings showing that hGR variants 1 B , 
1 C and 1 H are highly expressed in peripheral blood whereas hGR 
1 F is expressed almost exclusively in the brain, 23 and suggests that 
modifications in peripheral hGR 1 B and 1 C expression may be 
associated with the pathophysiology of PTSD. It is still unclear 
whether these changes in hGR expression reported in the blood 
might also be found in the brain. For instance, similar transcrip- 
tional differences in the hippocampus (HPC) may alter the 
inhibitory feedback exerted on the HPA axis activity that would 
consequently associate with the hypoactivity of the axis reported 
in this study. More work is required to define whether changes 
occurring in the peripheral tissue may also be found in the brain. 

The expression of specific hGR non-coding exon 1 variants has 
been associated with different GR protein isoforms. 49 These 
proteins include GRa, GRp and GR-P. GRa acts as the main 
functionally active isoform while GRP and GR-P inhibits and 
enhances GRa and p activity. 50-54 Interestingly, GRa and GR-P 
have been associated with hGR 1 C and 1 B expression, 
respectively 49 Thus, our results suggest that both GRa and GR-P 
may be up-regulated in individuals with lifetime PTSD providing a 
potential explanation for the hypoactivity of the HPA axis found in 
these patients. In addition, we found a correlation between hGR 1 B 
expression and Cortisol levels at awakening, 30 and 45 min after 
awakening. This suggests that hGR 1 B variant may be involved 
more specifically in the morning increase in Cortisol release. 
Indeed, one may assume that, under normal conditions, the 
enhancing effect of GR-P may be attenuated to allow more 
Cortisol to be released. However, the overexpression of hGR 1 B in 
PTSD, leading to higher levels of GR-P, may counter this 
attenuation by potentiating the GR-induced inhibition of the 
HPA-axis activity and subsequently decreasing morning Cortisol 
release. More work is required to clarify the roles of GR protein 
isoforms on Cortisol release. 

Our findings suggest that hGR expression may be regulated by 
changes in DNA methylation. DNA methylation is an epigenetic 
mark associated with transcriptional repression that acts by 
interfering with the binding of transcription factors and with the 
recruitment of the transcriptional machinery to the DNA. 16 Our 
results show a significant overall and site-specific hypomethyla- 
tion in the promoters of hGR 1 B and 1 C in the PTSD group 
compared to controls. In addition, hGR total and hGR 1 B 
expression was inversely correlated with overall and site-specific 
methylation levels in hGR 1 B promoter. This suggest that lower 



levels of methylation in hGR 1 B promoter may allow the 
transcriptional machinery to bind within this region and 
potentiate expression as it has been suggested for multiple hGR 
variants' promoters. 55 

Previous studies suggested that epigenetic mechanisms may be 
involved in the pathophysiology of PTSD. For instance, DNA 
methylation alterations in the promoter of the brain-derived 
neurotrophic factor gene associated with changes in gene 
expression, have been found within the HPC in a rat model of 
PTSD. 56 In humans, Uddin et al. 57 used a genome-wide approach 
and reported associations between PTSD and methylation 
patterns in peripheral blood noting that differentially methylated 
CpG sites were overrepresented in genes related to immune 
function and inflammation. This was associated with overall DNA 
methylation levels being inversely correlated with gene expres- 
sion levels. Collectively, these results suggest that traumatic 
events may induce alterations in DNA methylation patterns, which 
may change gene expression patterns and subsequently induce 
behavioral modifications in PTSD. It should be mentioned that we 
did not measure GR protein levels in this study and, although 
mRNA levels are consistently correlated with DNA methylation 
levels, future studies should consider measuring protein levels to 
strengthen the functional impact of DNA methylation on 
biological functions. 

Early-life stress has been shown to interfere with gene 
expression by altering DNA methylation patterns at multiple 
genes including GR. 1 7-19,56,58-68 For instance, decreased GR 1 7 
expression associated with promoter hypermethylation has been 
reported in the HPC of low licking and grooming mothers' 
offspring. 17 Similarly, lower hGR 1 B , 1 C and 1 F expression has been 
associated with promoter hypermethylation in the HPC of abused 
suicide completers 18,19 while hGR 1 B , 1 c , Id/ h, 1h ana " h 
expression patterns have been shown to vary in the brain of 
depressed individuals with no differences in promoter methyla- 
tion levels 40 However, our results suggest that trauma exposure 
together with a lifetime diagnosis of PTSD is associated with 
hypomethylation of hGR 1 B and 1 C promoter. Although contrast- 
ing, these results are in agreement with the variations in the 
activity of the HPA axis found in victims of child abuse and 
depressed individuals compared with lifetime PTSDs. Indeed, the 
HPA axis in victims of child abuse, as well as depressed patients, is 
hyperactive 69 " 71 In contrast, our results, as well as previous 
findings, 13,15,41,42 show that the HPA axis in PTSD patients is 
hypoactive and matches with the overexpression of hGR 1 B and 1 C 
associated with promoter hypomethylation found in PTSDs. 
Consequently, this suggests that the consequences of stress on 
the epigenetic regulation of hGR may differ according to type and 
timing of stressful events and psychiatric condition. Unfortunately, 
the design of our study did not allow to adequately control for 
depression. As such, it is unclear whether or not depressive-like 
symptoms may partly account for our results. More work will be 
needed to shed light on the important contribution that depres- 
sion has on the molecular characterization of PTSD symptoms. 

The regulation of the HPA axis takes place mainly in the HPC via 
the binding of GR by glucocorticoids inducing an inhibitory 
feedback on the activity of the axis 9,71 However, given that the 
participants for this study are alive, our experiments were 
performed in peripheral blood samples, which have been 
frequently used in studies of PTSD. 14,48,57 Individual differences 
in lymphocyte GR responsiveness have been associated with 
changes in the HPA-axis activity 72 and with PTSD symptoms. 73 In 
addition, T lymphocytes have a major role in cell-mediated 
immunity and it is believed that traumatic events may induce long 
lasting epigenetic marks conferring this cell type a memory of past 
adverse life events. Consequently, we believe that these cells may 
represent a proxy to the molecular alterations that occur in the 
brain. More work is required to determine whether similar 
alterations can be observed in the blood and the brain. 
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It is important to mention that our analyses revealed no change 
between current and lifetime PTSD. This may be owing to the 
small sample size of this exploratory study. Alternately, it is 
possible that these epigenetic changes, which may be induced by 
traumatic events, are stable in time and may not be affected by 
treatment. These changes may also increase vulnerability to 
stress in individuals with a past PTSD. Epigenetic mechanisms and 
gene expression in the brain and in peripheral tissue can be 
affected by antidepressant treatment. 65,66,74 However, our results 
suggest that this may not be the case for the regulation of 
GR in PTSD. This is also consistent with our previous results 
obtained from the brain tissue of abused suicide completers. 19 
Thus, although we did not correct for medication in this 
study, our analysis comparing current and past PTSD suggest 
that medication may not have an effect on Cortisol levels, GR 
expression and DNA methylation. Nevertheless, more work is 
required to study the effect of treatment on epigenetic mech- 
anisms and on the epigenetic control of GR expression in 
the blood. 

To summarize, our results suggest that the expression of hGR 1 B 
and 1 c is increased in peripheral blood samples from individuals 
with lifetime PTSD. These differences in expression were 
associated with overall and site-specific DNA hypomethylation in 
T lymphocytes and with lower release of morning salivary Cortisol 
in PTSD. Thus, our results suggest that traumatic events in PTSD 
induce DNA methylation alterations that modify gene expression 
and the activity of the HPA axis. However, these alterations differ 
from the ones described in depressed and abused suicide 
completers suggesting that different epigenetic mechanisms are 
involved in the regulation of GR in PTSD individuals. More work is 
required to better understand the epigenetic mechanisms 
regulating the expression of GR and the activity of the HPA axis 
in PTSD. 
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